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1. INTRODUCTION 

In recent 15-20 years, a theoretical investigations of properties of quark-gluon plasma 
has been of great interest. It is connected with intensive looking for a QGP in the 
experiments with collision of ultrarelativistic heavy ions. 

Two methods to study of the nonequilibrium phenomena in a QGP are used: method of 
temperature Green functions and kinetic approach. Significant progress has been achieved 
in the development of the first method. The effective perturbative theory was constructed 
in the papers by Pisarski and Braaten, Frenkel and Taylor [1] on the basis of the resum- 
mation of so-called hard-thermal loops (HTL's), and the problem on the sign and gauge 
dependence of the damping rate of the long wavelength excitations of QGP was solved. 
The progress of the thermal QCD makes possible a new view at the existing of kinetic 
theory of QGP, the basis of which was laid in the papers [2] by Heinz, Elze, Vasak and 
Gyullasy, Mrowczynski and others, and it has given impetus to its further development. 

In spite of the fact that the language and methods appear very different, there are close 
similarities between HTL approach and transport theory. Originally the kinetic theory 
was used by Silin to derive HTL in the photon's self-energy [3]. HTL's in the quark and 
gluon self-energies can be computed similarly. Moreover, Kelly, Liu, Lucchesi and Manuel 
[4] have shown that the generating functional of HTL's (with an arbitrary number of soft 
external bosonic legs) can be derived from the classical kinetic theory of QGP. This points 
clearly to the classical nature of the hard thermal effects. 

A further step in development of kinetic theory was made by Blaizot and lansu [5]. In 
contrast to the early papers on transport theory of QGP [2], these authors use from outset 
the ideas developed in thermal QCD in deriving of the kinetic equations. The equations 
obtained by them isolate consistently the dominant terms in the coupling constant g in 
hierarchy of equations which describe the response of plasma to weak and slowly varying 
disturbances, and encompass all HTL's. However, here it should be noted that if the 
influence of the average fermionic field is neglected, then the expression for current induced 
by soft gauge fields, obtained in [6] (and the nonlinear equation of motion, connected with 
it) fully coincide with the corresponding expression obtained in [4] from usual classical 
kinetic theory on the basis of consistent expansion of distribution function in powers of the 
coupling constant. In spite of the fact that intermediate approximation schemes in which 
these equations were derived, mix leading and nonleading contributions with respect to 
the powers of g and so, are not entirely consistent, this, somewhat justifies use of the 
classical and semiclassical kinetic equation found in [2]. 

Such close interlacing of two methods of investigation of nonequilibrium phenomena 
in QGP leads to the question: can one calculate whether the damping rate of bosonic 
modes corresponding to the hard thermal result [7], remaining whithin the framework of 
classical (semiclassical) kinetic theory only? Xiaofei and Jiarong [8] were the first to put 
this question. Because of obtained results, they have given a positive answer. 

As was shown by Heinz and Siemens [9] , in linear approximation the Landau damping 



is absent in QGP. In fact, the only mechanism, with which one can associate the damping 
following from the kinetic theory with damping from HTL approach, is so-called non- 
linear Landau damping. It bound up with the nonlinear effects of interaction of waves 
and particles in QGP. The multiple time-scale method which has proved successful in 
studying the nonlinear properties of electromagnetic plasma (EMP) [10], was used in [8] 
for determination of this association. By means of this method the nonlinear shift of an 
eigenfrequency of longitudinal oscillations in the temporal gauge, the imaginary part of 
which defines required nonlinear Landau damping rate, had been obtained by Xiaofei and 
Jiarong. Futher, the limiting expression of the derived damping rate for k = 0-mode was 
obtained, and numerical computations for approximate estimate were performed. The 
value derived by this means is in close agreement with similar numerical one obtained by 
Braaten and Pisarski [7] on the basis of effective perturbative theory. 

However, under close examination of above-mentioned paper we found certain inac- 
curacies in computations, which were of both principle and nonprinciple character. The 
elimination of these inaccuracies finally leads not only to a numerical modification of 
the limiting value of nonlinear Landau damping rate obtained in [8], but what is more 
important, it changes the sign of obtained expression (this subject will be considered in 
detail in Sec. 10 below). This points to some prematurity of the statements in [8] on ob- 
tained connection between nonlinear Landau damping rate and damping rate from the 
HTL-approach. 

Moreover, the physical mechanisms lying in the basis of the process of nonlinear scat- 
tering of longitudinal waves by QGP particles has not been revealed by Xiafei and Jiarong. 
The contribution to this process from the effect of self-action of the gauge field was 
dropped, although in a " soft" region of excitations it is of the same order as the contri- 
butions, bound up with the effects of the medium. An important problem of the gauge 
dependence for the obtained expression of nonlinear Landau damping rate was not con- 
sidered. In connection with above remarks it is evident that this problem requires further 
investigation. 

In this paper, we consider the above problem, using the approach based on obtaining 
of generalized kinetic equation for waves in quark-gluon plasma developed by Kadomtsev, 
Silin, Tsytovich and others [11] in connection with EMP as a basic method of investigation 
of nonlinear processes in QGP. 

We have shown that within the limits of nonlinear theory, developed for EMP and im- 
mediately (i.e. without any crucial changes) applied to investigation of nonlinear processes 
in QGP, nonlinear Landau damping rate 7'(k) defines not an inverse time of damping at 
the expense of absorption (dissipation), but the inverse time of spectral pumping of the 
energy of waves in the direction of a small wave numbers. The inequality 7^(0) < is 
a consequence of this fact, i.e. k = 0-mode, in contrast to [8], in this approximation 
increases. 

The outline of the paper is as follows. In Sec. 2, we derive a system of self-consistent 



equations in the covariant gauge for regular and random parts of the distribution func- 
tions of QGP particles and gauge field is obtained. In Sec. 3 the first order approximation 
of the color current is considered and the correlation function of random oscillations is 
introduced. In Sees. 4 and 5 the second and the third orders approximation of the color 
current are studied, and the terms leading in the coupling constant are separated. In 
Sec. 6 the generalized kinetic equation for waves in QGP is derived. From there the ki- 
netic equation describing the nonlinear interaction of longitudinal oscillations is extracted 
in Sec. 7. In Sec. 8 the physical mechanisms defining the nonlinear scattering of waves by 
plasma particles are considered. In Sec. 9 the estimate of a value of the nonlinear Landau 
damping rate for longitudinal eigenwaves is made. In Sec. 10 the inaccuracies, which were 
made in [8] in the process of computation of nonlinear Landau damping rate in the tem- 
poral gauge, are indicated, and its correct expression is written out. Comparison of this 
expression with similar one, obtained in the covariant gauge is performed. In Conclusion 
possible ways of further development of the scrutinized theory are discussed. 

2. THE INITIAL EQUATIONS. THE METHOD OF AVERAGING OVER 
STATISTICAL ENSEMBLE 

We use metric g'^'^ = diag{l, —1, —1, —1) and choose units such that c = ks = 1- The 
gauge field potentials are Nc x A^c-matrices in a color space defined by A^ = A"t" with 
N^ — 1 hermitian generators of SU{Nc) group in the fundamental representation. The 
field strength tensor Ff^^, = F^j^t"" with 

F;, = d,Al - d^Al + gr'^A'^Al (1) 

obeys the Yang-Mills (YM) equation in a covariant gauge 

d.F^^'^ix) - ig[A^{x), F^^{x)] - C'd^d^A^{x) = -fix), (2) 

where ^ is a gauge parameter, j" is the color current 

f = ge J d'pp'^lSpt^if^ - /,) + Tr (TV,)], (3) 

where T" are hermitian generators oi SU{Nc) in adjoint representation {(T"')'"^ = —if"''^'^, Tr(T"T^) 
NcS"''^). We denote the trace over color indices in adjoint representation as Tr. Distribu- 
tion functions of quarks fg, antiquarks fq, and gluons fg satisfy the semiclassical kinetic 
equations (neglecting spin effects) 

1 df - 

P'1^,U,±^9P''{F,.,^} = 0, (4) 

1 df 



where V^ and V^ are covariant derivates which act as 

V^ = df,-ig[Af,{x),-], 

V,, = d^-ig[Af,{x),-], 

[ , ] denotes commutator, {, } denotes the anticommutator, and A^, J-'^u are defined as 
A^ = A'l^T"', Tp,u = F^^T"". Upper sign in the first equation (Q) refers to quarks and lower 
one - to antiquarks. 

We begin with consideration of kinetic equations @. The distribution functions fq^q 
and fg can be decomposed into two parts: regular and random (turbulent) ones 

fs = f^ + fl . s = q,q,g, (5) 

so that 

ifs) = /f , ifl) = 0, (6) 

where angular brackets (■) indicate a statistical ensemble of averaging. The initial values of 
parameters which characterize the collective degree of a plasma freedom is such statistical 
ensemble. For almost linear collective motion to be considered below it may be initial 
values of oscillation phases. 

Further we set 

A, = A^ + Al, {AD = 0, (7) 

by definition. For simplicity the regular part of the field Aj^ will be considered equal to 
zero. 

Averaging the equation (^ over statistical ensemble, in view of (i|)-(0), we obtain the 
equations for the regular part of the distribution functions f^^ and f^ 

T 1 PifR 



1 9r 1 df^ 



1 dF 

T-^9P'{{{FL)nl,4z:^]), (8) 



p'dX = ^9P'{[Al, O - ::9P'{{{:fI)l, -^}) - -,9P' {{{:fI) n l) . -^}- 



1 df^ 1 df^ 



Here, indices "L" and "NL" denote the linear and nonlinear parts with respect to field 
A^ of the strength tensor (P . 

Subtracting (H) from (^), we define the equations for /^^ and fj 

P'd.flq = zgp^iiAD fl^] - {[AD Z^)) T 19P'{{FD)l, ^}t 



1 df^ df^ 1 df^ 



l9P'i{{F^.)NL,^} - {{iF^.)NL,^ 



1 dF dF 

T -^gp'{{{FL)NL, -i^} - {{{FL)nl, -i^])\ (9) 



1 df^ 

Ig — "yf Kl-^ii^Jg i- M-^ij.^ Jg i/J - 29P^{('Ffiu)L,-Q- 



P'd.f = ^9P^{[A'f] - {[A'f])) - T^gP^i^pL,^} 



1 dF dF 1 df^ 

—^gp^im^V, -^} - ({(^i)L, ^})) - ^gp^CFlUL - (Cf^nl), -^}- 

1 df^ df^ 



The system of equations (P) and (P) is suitable for investigation of nonequilibriuni pro- 
cesses in QGP such that the excitation energy of waves is small quantity in relation to 
the total energy of particles. In this case it can be used expansion in series in power of 
oscillations amplitude of the random functions /J 

oo 

/J = E/J^"^^=5'5"'^' (10) 

n=l 

where index n depicts that /J*-"^ is proportional to the n-th power of Aj^. Substituting 
expansions ( p!0D into @, and collecting terms of the same order with respect to Aj^, we 
derive the system of equations 

p'dX^'' = Tlgp'{{FpL,^h (11) 

p'dJ^P = ^gp''{[Al, fJ,P] - {[Al /J«])) T (12) 

1 df^^^^ df^^^^ 1 df^ 



p'dXP = ^gp'i[Al, /Jf ] - {[Al /Jf ])) T i^P^({(Fl)L, %^} - (13) 



-{{{FDl, ^})) T V({(^I)iVL, ^} - ({(Fjj^., ^})) etc.. 



Similar equations are obtained for /J'-"^ , n = 1,2,3, 
The nonlinear color current is expressed as 



oo 



J/J- Jfi 'J/J. ' WM/ J/j. 5 J/^ Z^Jf^ ' l-'-^/ 



n=l 



where 



jj(") = gf^ I d'pp.lSpf^ifJ^-^ - /J(")) + Tr (TVJ("))]. (15) 



Now we turn to the Yang-Mills equation (Q), connecting the gauge field with the color 
current. Averaging Eq. (§) and subtracting the averaged equation from (|]) in view of 
Eqs. (^ and (|T|), we find (for A^ = 0) 



d^iF^^L - r'd'd^Al + f^'> = -3l\ + igd^{\A^^, A^^\ - ([A^^ A^^)) + (16) 

+^^([^L {F^nL\ - ([Aj, {F^nL\)) + 9\\Ai^ [A^^ a^i] - ([^j, [a^^ a^^^ 



Here, in the left-hand side we collect all linear terms with respect to A and we denote: 



M 



2^1^ = j^(2)i/_|_jT(3)i'_j_^ ... To account for nonlinear interaction between waves and particles 
in QGP (in first non- vanishing approximation over the energy of waves), it is sufficiently 
to restrict the consideration to the cubic terms with respect to oscillation amplitude in 
expansion (plQl). 



We introduce the following assumption. Eqs. (^ represent the transport equations for 
averaged distribution functions. The correlation functions in the right-hand side of these 
equations have meaning of the collision terms of QGP particles with waves and describe 
the influence of plasma waves to a background state. 

We suppose that a characteristic time of nonlinear relaxation of the oscillations is 
small quantity as compared with a time of relaxation of the distribution particles f^. 
Therefore we neglect by variation of regular part of the distribution functions with space 
and time, assuming that these functions are specified and describe the global equilibrium 
in QGP 

fR ^fo ^ 2^IMMs(r,^) I f^ = fo = 2^^^6(v') ^ (17) 

Jq,q-Jq,q ^ (27r)3 ^-^ ^ e((P")T^)/^ + 1 ^ ^ ~ ^ (27r)3 ^-^ ^ e(p«)/^ - 1 ' ^ ' 

where Nf - being the number of fiavours for massless quarks, m^ is the four-velocity of 
the plasma at temperature T, and yU is the quark chemical potential. 

3. THE LINEAR APPROXIMATION. THE CORRELATION FUNC- 
TION OF THE RANDOM OSCILLATIONS 



Let us turn our attention to obtaining of kinetic equation for waves. The initial 
equation is Eq. (K). The left-hand side of Eq. (^) contains a linear approximation of 



the color current, explicit form of which is easily defined from Eq. ([Tl|). We prefere to 
work in momentum space; the corresponding equations are obtained by using 



A^{x)= d^kA^{k)e 



ikx 

1 



and similar translations for fj^, fT. The result of Fourier transformation for Eq. ( ]TT| ) is 
fq,q {k,P) = T9^-^. TT^A^ik) , f ^ >{k,p) = -g^-^. — —\Au{k), (18) 



Here, x'^^{k,p) = {pk)g''^ —p'^k^, and the suffix "T" for a gauge-field is omitted. Sub- 
stituting ([T8|) into (|15D (more precisely, in Fourier transformation of ([T5|)) we define a 
well-known form [2, 12] of linear over field approximation of the current 

/(i)'^(A;) = W{k)A,{k), (19) 

where 

p^'ip'^ikdp) - {kp)d;)Me, 



W{k) =g^ f d^p- 



'eq 



pk + ipoe 



is the high temperature polarization tensor, and J\feq = |(/° + /?) + N^fg. 

Further we rewrite Eq. (|IBp in the momentum space. Taking into account (^^j, we 
obtain 

[k'g'^^ - (1 + r')k^k'' - W^ik)]Alik) = jltik) + 

+ f'^' J sifXi^tikMiih) - {Al{k,)Ai{k^)))5{k -k,- k2)dk,dk2 + (20) 

+/'=^/^'^/sr.':..,.3(^e(^i)^i(^2)A:(A:3) - (A:(fco^i(^2)A:(A:3))) 

S{k — ki — k2 — k^) dkidk2dk3, 
where 

sitt = -^9{k^9'' + k^29''' - ki^g""-) , ^^tMM = d^g^'g'"- (21) 

Let us multiply Eq. (|^) by the complex conjugate amplitude A*^{k') and average it 
[fcV^ - (1 + r')k>^k'' - Yi^^{k)]{Al\k')Al{k)) = {A*;{k')jlt{k)) + 

+ f^ J SifX{A*;ik')Kiki)Ai{k2))6{k -k,- k,)dk,dk, + (22) 

+f'ff^'' j^';,tMM^K'^k')^'i^^kMi{kMl{k,))5{^ -k,-k,- k,)dk,dk,dks. 

We introduce the correlation function of the random oscillations 

i;t{k',k) = {A*;{k')Al{k)). (23) 

In conditions of the stationary and homogeneous of QGP, i.e. when the correlation func- 
tion (|23|) in the coordinate representation depends on the difference of coordinates and 
time Ax = x' — x only, we have 

I^t{k',k) = i;t{k')6ik'-k). (24) 



By the effects of the nonhnear interaction of waves and particles, the state of QGP 
becomes weakly inhomogeneous and weakly nonstationary. The dependence of I°jl on 
k' — k = Ak is no longer ^-shaped, i.e. k' ^ k and it "smeared" on Ak, with | Ak/k \<^ 1. 

Let us introduce I^l{k',k) = I^^{k, Ak), and insert the correlation function in the 
form 

i;l{k,x) = Ji;l{k,Ak)e-''''^dAk, 

slowly depending on x. In Eq. (521) we make the change k ^ k' , a ^ b, complex 



conjugate and subtract obtained equation from Eq. (p2D , beforehand expanding of the 



polarization tenzor into Hermitian and ant i- Hermit ian parts 

W'^ik) = n^^"(A;) + n^^"(A;) , n^'^"(A;) = U*^'"'{k) , n^'^"(A;) = -U*^'"'{k). 

By using the definition (|23|), we obtain the result 
[{eg^"" - (1 + r^)k^'k'') - {k'^g"'' - (1 + C^)k'^'k"') - (n^^^(A;) - U^'"'{k'))]lf^{k' , k)- 

-[Ii^^\k) + Ii'^^'^{k')]I';^^{k\k) = {A*;{k')jlt{k)) - {Al{k)f^rik')) + 
+1'"^ J SifX{A*;{k')At{k^)Ai{k,))6{k -k,- k,)dhdk,- 

- r' J SliZl'{4{k)A:^{k^)At{k,))6{k' -k,- k,)dk,dk, + (25) 

+f'^ f^'"" j ntMM^K'^k')^l^kMi{k2)Al{k^))5{k -k,-k2- ks)dk,dk2dk,- 

-r^^/^'^/srfcJr.„.3(4(^)^r(^i)^i'(^2)Ar(A;3))5(A;' -k,-k,- k,)dk,dk,dk,. 

We assume that anti-Hermitian part of 11"^ is small in comparison with 11^ and it is 
a value of the same smallness order, as the nonlinear terms in the right-hand side of Eq. 
(p5|). Therefore it can be suggested that Ii'^'^'^{k) ^ n^^°'(A;'), and the term with 11"^ can 
be rearranged to the right-hand side of Eq. (|25|) .The remaining terms in the left-hand side 
of Eq. ( p5D we expanded in a series in powers of A A; to first smallness order. Multiplying 
obtained equation by e"*^'^^ and integrating over A A; with regard to 



f Akxi;t{k,Ak)e-'^''^dAk 



dl^t{k,x) 
i- 



we obtain finally 

[k'^g^"' - (1 + c^)k''k'' - u^'^'ik)]—^ = 2m^'"'rl- 



dkx' V ' s 7 ' ^ '^ dx^ """ ^^'^ 

-z fdk'{{A*;{k')jlt{k)) - {A'^{k)fJr{k'))} 



- r'j dk'dhdk, siyi^;'^{A'^ik)A:%kMtihmk' -k,- k,)} - (26) 

-^{/'"^^/^'^/s^:fct..,.3(^r(^')^'.(fci)^l(fc2)A:(A:3))5(A: -k,-k,- k,)dk,dk,dk,- 

-rf ff'^ j n''t^.MK^k)^i''^kMi\kMi\h))5{k' -k,-k2- k,)dk,dk2dk,}, 

where jlTik) = jr(2)«/^(A;) + j^(3)aM(^), 

We made several remarks relative to obtained Eq. (]26|). The term with 11^ introducing 
in the right-hand side of Eq. (|26|) corresponds to linear Landau damping. However, as 
was shown by Heinz and Siemens [9] , linear Landau damping for waves in QGP is absent 
and hence this term vanishes. 

As it will be shown below, the terms with E make no contribute to processes of the 
nonlinear interaction of waves with particles of QGP, and hereafter they will be dropped. 

4. THE SECOND APPROXIMATION OF THE COLOR CURRENT 



Now we conserned with computation of the nonlinear corrections to the current in the 
right-hand side of basic Eq. (^). In this section we consider the second order oi j^^'^^^{k). 

At first we define /g J . We carry out the Fourier transformation of Eq. (0) 
-^{pk)fl^'\k,p)=gp^J {z{[A^{k^)jl^'\k2,p)] - {[A^{k,)jJ}'\k„p)]))T 

±^gi[A^ik^),A,{k2)] - {[A^ik^), A,{k2)]))^}Sik - fci - k2)dk,dk2. 

Substituting the obtained fg^^ from (|18|) into the last equation and collecting similar 
terms, we obtain 



.T(2)_ .[^^i1pV / (^2^yjf),,6,,. ^,c,,. ^ fAb, 



f,r = Tg' ' 12 / \ ''''" {Al{k,)Al{k,)-{Al{k{)Al{k,)))5{k-k^-k,)dk^dk,+ 

pk + tpoe J pk2 + ipQt 



9' {t\t'} f..uxn. ^^ 9 ^X'^^(fc2,pX,w.^^.^.c„^ , ,b 



(27) 



i. ^ • / x'\h,p)Tr^ \ Y'^' -^ {Al{k,)Al{k,) - {AlikMUh))) 
2 pk + ipQt J op^ \pk2 + ipoe opp j " '^ 

5{k — ki — /C2) dkidk2. 



10 



The expression for fT^'^\k,p) is obtained from (^) by choosing upper sign and replace- 
ments fg — > fg, t"" — >• T". Substituting obtained expressions /J*-^-*, s = q,q,g into (JT^) (for 
n = 2), we find required current correction 



7 pfc + ipQt pk2 + tpoe 

6{k — ki — k2) dkidk2+ 

4^ I ,',p^^m^A (P^'-^^^) iAiikM%k.)~iAiikMiim) 



4 J pk + ipoe dp^ \pk2 + ipoe dpP ^ 

6{k — ki — k2) dkidk2- 

The contribution of gluons to the expression with symmetric structure constant d"'^'^ here 
drops out. This is connected with the fact that in calculation of trace of anti-commutators 
we have: Spt'^jt'', t'^} = \d°''^^ -for quarks and antiquarks, and Tr T^jT^jT*^} = - for 
gluons. The symmetry of contributions can be restored if we note that besides usual 
gluon current ig{x) = gt"" J d^pp'^Tr{T°'fg{x,p)), the kinetic equation for gluons admits 
a covariant conserving quantity 

Xgt-J d'pTT{V^fg{x,p)), (29) 

where (V°-y"^ = d"'^'^ and A is a certain arbitrary constant. The covariant continuity of 
(^) is evident from the identity: [V.T^] = if'^'V^. On addition of (^ to (0) we have 
contributions to the nonlinear current corrections only. Adding (pOJ) to the second current 
iteration (|I5|) and taking into account the equality 



instead of (|28|) we derive more general expression for j^^^^ 

^■T(2)a,(^) = / 5,tr.'(4(fcl)^A(fc2) - {AlikMlik2)mk - k, - k2)dk,dk2, (30) 



wnere o^mm ~ J "^kMM + " "^kMM ' 



q(II)^.uX _ .3 [ M^ P^P'P^ {k2dpKg) .„, ^ 

^(iii),ux _9'' /■ ,4 P^X'^Jkup) d f x''{k2,p) dMeA ._. 



2 J pk + ipoe dp" \pk2 + ipoe dpP 

1 



-Me.= -(/°-/°) + AiVjO. 



The tensor structure of Sf, ^^ j^^ exactly coincides with appropriate expression ob- 
tained in calculation of j-^*^^)^ in electromagnetic plasmas [11], and hence part of current 

11 



with d""^^ has a meaning of Abehan part of the color current j^(2)a/i_ rjnj,^^ term with S*^ ^^'^^ 
is purely non-Abelian, i.e. it has no Abelian counterpart. 

Let us estimate orders of S*^^^-* and S^^^^\ Following usual terminology [1], we call an 
energy or a momentum "soft" when it is of order gT, and "hard" when it is of order T. 
We will be considered, as in Ref. [5], that collective excitations carrying soft momenta, 
i.e. k ~ gT, and plasma particles have the typical hard energies: p ^ T. On this basis, 
we have the following estimate for S^^^^ 



Q{II)fJ.UX 



g'T. 



(33) 



Here, we considered that by virtue of the definitions ([T 



A4g~ 



1/T2. 



In expression (|32|) the integral of energy with gluon distribution function is infrared 
divergent. In a similar manner [8], we regulate it by introducing an electric mass cut-off 
of order gT, and only take the leading term in g. Then it can be found that in (^), the 
part related to the gluon distribution function is of order g^T{g In g) and the other part 
related to the quark and antiquark distribution functions is of order g'^T, i.e. 



^iin)fiux 

^k,k-iM 



g'T{g\ng)+g'T. 



Hence 5*^^^^ which is purely non-Abelian, is of lower order in the coupling constant 
than S^^^^\ which has an Abelian counterpart. This fact was first seen in Ref. [8]. 

5. THE THIRD APPROXIMATION OF THE COLOR CURRENT 



Now we calculate the nonlinear correction of j"^*-^-''^. We carry out the Fourier transfor- 
mation of Eq. (|l^) and taking into account (|18D, (^), after cumbersome transformations, 
we define 



fl?\k,p) 



rhcf icde^f 



pi^p^pO 



ik2dXc 



-J- i J J f f t' V ^-pjq^q, 

pk + ipoe p{ki + k2) + ip^t pk2 + ipo 



tg^t'',{t'',t'}] 



p''X^''{h,p) d ^ X'"'{k2,p) dfl 
2 pk + ipQt p{ki + ^2) -|- ip^e dp°' ^pk2 + ipoe dpp 



q,Q 



+ 



+ — ^— X {h,P)- 



T 



2 pk + ipoe 
g'{tMt^,t^}} 



pXpa 



ihdX 



q,q^ 



dp'^ ^p{ki + ^2) + ipo^ pk2 + ipo'^ 



T 



d 



^^^(^-P)^( 



X^"(A;i,p) d fX'''{h.p)dfl 



{-, 



'Q,q 



4 pk + ip^e '^ ^ •^'■f ^ QpT \p(^ki + /C2) + ipo^ dp°' ^pk2 + ipo€ dpP 
{Al{k,)Aiiki)A:{k2) - Al{k,){A'i{kMlik2)) - {Al{kMiiki)A:{k2))) 



(34) 
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6{k — ki — k2 — ks) dkidk2dk^+ 

2 J pk + ipQt op'^ ^pks + tpoe opP ' 

6{k — ki — k2 — k^) dkidk2dk'i. 

The gluon part of fT^^\k,p) is obtained from ( ^4|) by replacements /° — > fg and t° — > T". 
Substituting obtained expressions for f^q and /J*^^-* into (|1^) (with regard to additional 



current (^Of )) and taking into account the equalities 

Tr ({T^ T^}{T^ T'^}) = N^'^^^d"'^' + 45"'',5"'^ + 25"'^^'^'' + 28''''8^\ (35) 

Tr({P",T^}{T'^,T"}) = 0, 
we find the required form of j^(3)«/^ 

^■T(3)a.(^) = I sSr,C/r3(4(A:3)Al(A:i)A:(fc2) -4(A:3)(Al(fci)A:(A:2))- 

- {Al{k^)Ai{kr)Al{k2)))5{k -ki~k2- k-i) dkidk2dk^ + (36) 

^^atcfcej RP^^\^,^^^^,^(Al{ks)AiikM:ik2) - {A^ik,) AHh) Al{k2))) 
6{k — ki — k2 — /cs) dkidk2dk^. 
Here, 

spabde^ivXa rabc fcde-sp{I)fJ.i^>'Cr cabcjcde-^(II)pi.v>^cF .i-abi:de-y{III)fJ.!^>^cr,mbcfcde-sp{IV)fiu\a 

^k,kiMM ~ J J ^kMMM "I" / " ^k,kiMM "•" " " ^k,kiMM "T " / ^A:,fci,A 



.(I)f.uXa _ 4 /■ .4 P^P'^P^P'' I {k2dpMeq) 



,k2,k3 



^(i)fiUAa _ 4 I J4^£^_r_r_r - yz^-p-' - eqj /^yN 

' ' ' J pk + ipoe p{ki + k2) + ipoe pk2 + ipoe 

, 1, 2, 3 2A^^ 7 pk + ipoe dp'^ ^p{ki + k2) + ipoc <9p" ^p/i;2 + ^Poc <9p'' ^ 

(38) 

y(V)iiu\a _ Ncy(III)fj.uXrT 
k,ki,k2,k3 r) k,ki,k2,kz ' 

The expression for S^^^^ is obtained from (^) by exception of quark and antiquark 
contributions. The availability of the term with S*^^^^ is reflection of more complicated 
color structure of the gluon kinetic equation in comparison with quark and antiquark 
equations. 

The terms with S*^^^\ S'^'^^^ and R are defined as the interference of Abelian and 
non-Abelian contributions. In the case of isotropic, homogeneous and colorless plasma, 
the correlation function ( pBj) is proportional to unit matrix in color space, i.e. plasma 
oscillations are degenerate in color association. This leads to the fact that the coefficients 
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standing before these interference terms, supporting the color indicies, are vanish in kinetic 
equation for waves and therefore their exphcit form is not given here. With the availabihty, 
e.g., of external color field only, when degeneration is removed (the correlation function 
(^) becomes nontrivial matrix in color space), the interference of Abelian and non- 
Abelian contributions will be presented. 



4tt:ks - /• (39) 



At the end of this section we estimate the order of S'^'^^ and T,^^^^\ It follows from the 
expression ( |57D that 

-'k,ki,k2,k:j 

Cutting off, as in the previous section, integration limit for the gluon distribution function, 
we find 

By this means, purely non- Abelian contribution of S*^^^ is of lower order in the coupling 
constant than Abehan - S^^^^), E^^) and E^^^l 

6. THE GENERALIZED KINETIC EQUATION FOR WAVES 

Let us consider the initial equation for waves (|26|) . Substituting obtained nonlinear 
corrections of the induced current by field PD| ) and ( P^D into this equation, and considering 
the terms of leading order in g only, we obtain 

-[k^g^" - (1 + C^)kf'k'' - U^^'^ik)]- "" - 



-2 



dkx dx^ 

^-ij dk'{f'''S^f^^,^{A*;{k')A:{k,)Al{k2))dk^dk26{k -k,~k2 

-r'sir,'^^,^{Ai{k)A:%kMt{k2))dk,dk,6{k' -k,- k,)+ 
+f''^f^''^^kl':ktkM7(k')At^^^^ 

dkidk2dks5{k - ki - k2 - k^) - (40) 

-r^f^'^K'^'tk'^kM.i^^^^^^^ 

dkidk2dk'i5{k' — fci — ^2 — A^a)}. 

TT„„„ c/^^^ — cWf^'^^ I QiII)^^vX 
neie, o^ ^.^ ^.^ — '-'k^M "•" '-'kMM ■ 

By virtue of weak nonlinearity, oscillations phases of a field weakly correlate among 
themselves. Therefore mean value of four random quantities can be approximately divided 
into product of mean values of two fields. For mean value of three fields this decomposition 
vanishes, and it should be considered a weak correlation of fields. For this purpose we use 
the nonlinear equation of a field (^), taking into account in the right-hand side of (|20| ) 
the terms of the second order in A only 

[A; V - (1 + r^)k^'k'' - W{k)]A''^{k) = 
14 



= r'^J SZ',kMlih)Al{h) - {Al{h)Al{h)))6{k ~h- h)dhdh. (41) 
The approximate solution of this equation is in the form 



- {Af^\ki)Af\k2)))5{k -ki- k2)dk,dk2, (42) 

where A|J')"(A;) is a solution of homogeneous Eq. ( |4T| ) corresponding noninteracting fields, 
and 

V^^ik) = -[Pg^, - (1 + r')Kk. - ^^.AkT' (43) 

represents the medium modified (retarded) gluon propagator. 



Now we substitute ( ^2] ) into correlators of three fields introducing in Eq. (|4D|). Because 
^(0) represents amplitudes of entirely uncorrelated waves, the correlation function with 
three A^^^ drops out. In this case it should be defined more exactly each of the terms 
in {A*^''{k')Al{ki)A'i{k2)) and {A'l^{k)Al%ki)Af{k2)). In the correlation functions of four 
amplitudes, within the accepted accuracy, it can be done distinction between the fields A 
and /l(o). 



Finally Eq. (|40| ) becomes 



a grab 

— [k^g^" - (1 + r^)k^'k'' - U^^^'ik)]^^ 



dkx dx^ 

= -zj dk'dk,dk2dk,{f'^ff^^^5{k-k,-k2-k,)t^,f^^,^^,^{{Al^{^^^ 

-{A*;{k')A:ik,)){Ai{kMUk2)))- 
-rfff^^Sik' -k,-k2- k,)tl^l'll^^,^{{A'^{k)Al%k,)Al\k,)A*:{k2))- 

- Kik)AT{k,)){Al\k{)A7{k2m + (44) 

+if'Tf j dk' j dk^dk2dk[dk',{v;^{k')-v^,{k))s^%i,^si?tr^,,^^^^^ 

A*{{k[)A*J{k',)) - {Alik,)Aiik2)){A*{ik[)A*Jik',)))6ik - k, - k2)6{k' - k[ - k',). 
Here, 

Y^/iu\a _ ^piI)iiuXa- _ j Qjiup _ qfipi/ \^ ( U J- h \ caXa (Ar\ 

^kMMM ~ k,kiMM Wfc,fc3,fci+fc2 '~'kM+k2,k3> P'^^^^ '^> ki+k2MM' ^ > 

It follows from the definition (|43|) that the propagator is of order ~ \j g^T"^ . Taking 
into account (^) and (P^ ) we see that all terms in the right-hand side of ( ^ij) are of the 
same order. This explains why in the expansion of the current ([M]) the following term - 
jT(3) gj^ould be retained in addition to the first nonlinear correction j^^^-*: it leads to the 
effects of the same quantity order. 
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7. THE KINETIC EQUATION FOR LONGITUDINAL WAVES 



Let us divide the mean of four fields in Eq. ( PD into three pairwise products of 
two-point correlations by the rule 

{A{kM{k2)A{h)A{k,)) = {A{h)A{h)){A{h)A{h)) + {A{h)A{h)){A{h)A{h))+ 

+ {Aih)Aih)){Aik2)Aih)). 



Taking into account that the spectral densities in the right-hand side of Eq. ( PD can 
be considered as stationary and homogeneous those, i.e. having the form (P^, and setting 
jab ^ ^abj^^^ ^g gj^^ instead of Eq. (||) 

8 8T 

— [fcV^ - (1 + r')k'k^ - ^""-"m^ = 

= 2iV, J dk, Im{f:^fl^_,^ - f:tZA-JWk)IuAki) + (46) 

+NMiv,o.{k)) I dk.dk, {s',%1,^ - six^,;){sifx - s:jXy,x{k,)iUk,)Sik -k,- k,). 



As it is known [13, 14], in global equilibrium QGP the oscillations of three typies can 
be extended: the longitudinal, transverse and nonphysical oscillations. In this connection 
we define the spectral density I^u{k, x) = I^^ in the form of expansion 

V = P,A + Q,A + iD^Jl , 4*'''^^) = /(*'''") (fc, x), (47) 

where the transverse and the longitudinal projectors [14] are P^i, = g^y — k^k^/k'^ — 
Qfiu , Qfiu = u^Uy/u^, and u^ = k'^u^ - kf,{ku)] D^^ = kf.k^/k'^. By using these projectors 
the propagator (^) can be written in more convenient form 



Here, 11* = ^U^^'^P^^, 11' = H'^'^Q^i/- At finite temperature, the velocity of plasma intro- 
duces a preferred direction in space-time which breaks manifest Lorentz invariance. Let 
us assume that we are in the rest frame of the heat bath, so that m^ = (1, 0, 0, 0). 



Eq. (^Bp with the expansions ( ^Tj ) and (^) enables us to investigate various nonlin- 
ear processes in QGP: the nonlinear scattering of longitudinal waves in longitudinal or 
transverse waves; the scattering of transverse waves in longitudinal or transverse waves; 
the merger of two longitudinal waves in one transverse wave etc.. In this work we restrict 
our consideration to investigation of most simple process - the nonlinear scattering of 
longitudinal waves by particles of QGP in longitudinal those. 
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Further derivation of kinetic equation for longitudinal oscillations is similar to cor- 
responding derivation in the theory of electromagnetic plasma, therefore we restrict our 
consideration to its schematic description. 

Now we omit nonlinear terms and anti-Hermitian part of the polarization tensor in 
Eq. (P^D. Further substituting the function 5°-^Q^y{k)ll5{k' — k) instead of I^l{k', k), we 
lead to the equation 

Re{e\k))li = 0. 

Here, we use relation: l—Il\k)/k'^ = e\k). The solution of this equation has the structure 

li = Il5iuj - 4) + I-^5iiu + 4) , 4 > 0, (49) 

where /^ is a certain function of a wave vector k and cj^ = ^^'(k) is a frequency of the 
longitudinal eigenwaves in QGP. 

The equation describing the variation of spectral density of longitudinal oscillations is 
obtained from Eq. (|i6| ) by replacement: /^^ —>■ Q^u{k)l\^, where /^ is defined by (^91). 5- 
functions in (^) enable us to remove integration over frequency and thus we have instead 
ofEq. (H) 

P Qj^ I dx^~ J ^ kiUin[l^fc,fc,fci-fci ^fc,fci,fc-fciJ + 

k 

POO r 

+Nc I dto dkid'k2li^-^I]^.^{Gk,ki,k2 + Gk-kiM +Gk,ki,-k-2 + Gk,-ki,-k2)uji=ujl ,0^2=0;' 5 
Jo J 1 2 

where 

GkMM = Ml^,a{k)){S',ZM - SiZM)^SlfX - SlJ^\;)Q,y^{k,)Q,^(k2)8(k -k,- k,). 

Let us consider the terms entering in the right-hand side of Eq. (^Up. The integral 
with the function Gk^kiM ^^ different from zero if the conservation laws are obeyed 

k = ki + k2, (51) 

These conservation laws describe a decay of one longitudinal wave in two longitudinal 
waves. However for a spectrum of the longitudinal oscillations in QGP, the equahties 
(pl|) do not hold simultaneously, no matter what the values of the wave vectors k, ki 



and k2 may be, i.e. this nonlinear process is forbidded. Therefore the integral with 
GkMM vanishes. Remaining integrals with G-functions differ from (^) in that some of 
the interacting waves are not radiated but absorped. They also vanish. 
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The expression 

V-'k,k,-kl,kl ~ ^fe,-fcl,fc,fel) ^Ml^JVl^O-lft^l) \uJ=UjI,UJI=UjI 5 (52) 

contains the factors 

l/{pk + ipoe) , l/{pki + ipoe) , l/{p{k - ki) + ipQe), 

by the definitions of functions entering in it. Imaginary parts of first two factors should 
be setting equal to zero, because they are connected with linear Landau damping of 
longitudinal waves (which is absent in QGP), and therefore the imaginary part of the 
expression (|52D properly introducing in Eq. (|50D will be defined as 

^ --5(4-<-v(k-kO). 



Im 



p{k - ki) + ipoe 



w=a;J^,a;i=(^j^ Po 



It follows that nonlinear term in the right-hand side of (^) with the function (|52D is 
different from zero if the conservation law is obeyed 

^lc-<-v(k-ki) = 0. 

This conservation law describes the process of scattering of longitudinal wave (plasmon) 
in longitudinal one by the particles in QGP. 

Let us consider in more detail the term in (|52D (see definition (|45|) ) with propagator 



T^pa{k — ki). By expansion (^Sp this propagator represents the nonlinear interaction of 
longitudinal waves with longitudinal ones through three types of intermediate oscillations: 
the transverse, longitudinal and nonphysical oscillations depending on a gauge parameter. 
The term with 

/ Ppa{k2) 



/Co 




w=a;J^,a;i=ajJ^ 



(hereafter ^2 = A; — A;i) in general, describes two fundamentally different nonlinear pro- 
cesses: 

1. if ^2 = (^k — ^ki) k — ki) is a solution of the dispersion equation k"^ — n*(fc2) = 0, 
then this term describes the process of merger of two longitudinal oscillations in 
transverse eigenwave; 

2. otherwise, it defines the process of nonlinear scattering of longitudinal waves in 
longitudinal those through the transverse virtual oscillation (for a virtual wave in 
distinction to the eigenwave, a frequency uj and a wave vector k are not connected 
with each other by the dispersion dependence: cj 7^ u;(k)). 

The equality k\ — Il\k2) = does not hold for longitudinal oscillations, as we see above 
and therefore, the term 




LU^UJ. ,101^0^1 



defines only the process of scattering of longitudinal waves in longitudinal those through 
the longitudinal virtual oscillation. 

Let us consider the contribution of nonphysical intermediate oscillations 



■^1 



By direct calculation one can show, that in contraction with tensor Dpa{k2) the complex 
factor l/{p{k — ki) + ipoe is reduced in the expressions with S'-functions. In particular it 
follows that contribution from the process of nonlinear scattering of longitudinal waves, 
connected with nonphysical intermediate oscillations drops out. Therefore the gauge 
parametor ^ is absent in the equation for longitudinal wave. 

The remaining terms with E are distinguished from above considered terms by a sign 
of fci, and describe the processes of simultaneous radiation or absorption by particles of 
two waves. The contribution of these processes is exponentially small in relation to the 
scattering process and therefore these terms are omitted. 

Summing the preceding and going from the function /^ to the function 

k 

having the physical meaning of spectral density of energy of longitudinal oscillations, we 
find from (^) the required kinetic equation for longitudinal waves in QGP 

^ + VL^ = -7{(^)}1^^, (53) 

where 

, _ dul _ ^^dRee\k)^ /^dRee\k). 
Vi, 



dk 



(9k / du 



is the group velocity of longitudinal oscillations and 



•^k ■" ^ki 



'-^fiuXa y-i/ii/crA 



k'^kj ^ duj '' ^ duji 



^^^k,k,-ki,ki ^fe,-fci,A:,fci)QMA(^)<5i/o-(«^ly 



iU=Uj'^,U)l=Uj!^ 



(54) 



presents the damping rate caused by nonlinear effects and being the linear functional of 
spectral density of energy. 

One can write (^) in the form which is more close to usual representation if the 
spectral density of number of longitudinal oscillations is entered 

Nl = Wi/ul 
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It fulfils the role of distribution function of a number of plasnions. Then instead of (|53| ) 
we have 

^ = ^ + V,^ = -7{iVjiV,. (55) 

8. THE PHYSICAL MECHANISM OF THE NONLINEAR SCATTER- 
ING OF WAVES 

Now we transform 7'(k) to the form allowing more neatly explain the physical meaning 
of the terms entering in the nonlinear damping rate. The first transformation of this type 
was proposed by Tsytovich for electromagnetic plasma [11]. 

By the definition S (^5[) we have 

i'^k,k,-ki,ki ~ ^fc,-fci,fc,fei)<5/iA(^)v!^(T(«;i) = {^k,k-ki,ki ~ ^A:,-A:i,fc,fci)<5MA(^)vi^(T(^l) + 
+ 7.2 _ HlfU ] y^kMM ~ ^k,k2,kjy^k2,k-ki ~ ^k2-ki,k>Qpa{k2)Qfix{k)Qua{kl). (56) 

Here, for simplicity in the last term in the right-hand side of Eq. ( ^61 ) the effects connected 
with existence of transverse intermediate oscillations are neglected. At first we consider 
the expression with S^^^. By virtue of definition ( P7D we have 

I™ (^fc,A;,-fci,fci "" ^fc,-A;i,fc,fci)vMA(^)v!^o-(^l) = 
~ k^kjk^kl J ^ P02 d{Pf^2)[k2C)pJ\eg)- 

The contribution to 7'(k) from this term may be introduce in the form 

- J w^{k, ki)iV^^po(fc2fi'pA4J ^Vki, (57) 



where 



<(k,ki) 



p 

Tl{lr.\\ -1 



2nN, fdRee\k)\-^ fdRee\ki)\-^ ,, i i „ , xm.s., i m2 .roN 

TkmVA^^) ..A^j^) ■ i(.L-.L-v(k-kO)|A=^(k.kOr. (58) 



UJi^U)\ 



A^(k.k0^^^^^i^^l^^»i^l^. (59) 

|k||ki| oj^. — (kv) 



To clear up the physical meaning of contribution (|57|) , it is convinient to compare 
it with appropriate contribution in the theory of electromagnetic plasma. In this case 
as shown in [11] this contribution describes the Thomson scattering of a wave uj\^ by 
particles: a wave uj\^ sets particles of plasma into oscillation and oscillating particles 
radiate a wave uj\^^. The corresponding function U7^(k, ki) presents the probability of 
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this scattering. As it was shown above in quark-gluon plasma for a soft long-wavelength 
excitations all Abelian contributions is at most glng times the non-Abelian ones and the 
basic scattering mechanism here, is essentially another (next we consider it briefly, the 
details will be published elsewhere). 

For revealing this mechanism we use the classical pattern of QGP description [2], in 
which the particles states are characterized besides coordinate and momentum by the 
color vector Q = (Q"), a = 1, . . . , A^^ ~ ^ ^-Iso. As was shown by Heinz [2] there is an 
intimate connection between the classical kinetic equations and semiclassical ones (^. 
Therefore in this case use of classical notions is justified. 

Let the field acting on a color particle in QGP represents a set of longitudinal plane 
waves 

Ki^) = J [Q,uik)Al%^^.^e^^-^-l'dk. (60) 

The particle motion in this wave field is described by the system of equations 

m^=9Q"F-|;^. (61) 

^ = -Sf'^^AlQ'. (62) 

Here, r is a proper time of a particle. The system (|61|), (|6^) is solved by the method 
of successive approximations - expansion in the field amplitude. A zeroth approximation 
describes uniform restlinear motion, and the next one - constrained charge oscillations in 
the field (0). With a knowledge of the motion law of a charge, the intensity of radiating by 
it longitudinal waves can be defined. In this case Eq. (^) defines the Abelian contribution 
to radiation, whereas (|62D - non-Abelian one and interference of these two contributions 
equals zero. The scattering probability calculated by this means, based on Eq. (B^) is 



coincident with obtained above (|58D . 

In this manner the contribution ( |5^ ) to 7'(k) are caused by not the spatially oscillations 
of a color particle, as it occurs in electromagnetic plasma, but the initiation of a precession 
of a color vector Q of a particle in field of a longitudinal wave (^) (Eq. ( p^ ) conserves 
the length of a color vector) . 



Let us consider now more complicated term in ( pq ) connected with S'-functions. By 
exact calculation, using the definitions (|2l| ) and (^Tj), it is not difficult to see that the 
following equality is obeyed 

isZik, - sZlkJ Mk)Mk2)Mh) = - {sCk,-k, - sz-k,,k) Uk2)u>.{k)uAh) ^ 

Then the contribution to 7'(k) from ^-functions can be represented as 

1 f^Ree^{k)^-l ^dRee^{ki] 



2/V f r]k N^ i ( dKee\k) yi / CKee'(, 



'^i='^Li 
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{k^klklY ye\k2) 



'>kM 



yK~KlK2)' 

Next we use the relation 



UJ^Lul, LJl^UJ. 



Taking into consideration the equahty 

2 

lme\k2) = -^ I d''pipof5ipk2)ik2d,K,) , 



one can write contribution from the first term in the right-hand side of (p^ ) to the nonhnear 
damping rate in the form 



w^{k, k^)Nl^po{k2dpKg) d'pdk,, (65) 



where 



w^p(k,ki) 



A^(k, ki) = A^'^)(k, ki) + A^*^^)(k, ki) = 



|k||ki|ki \k^ d{k2 



a;=a;; , a;i=a;; 



Here, A"^ represents the sum of two contributions. The first contribution defining by 
the function S^. l^ is connected with self-action of a gauge field and had no the analogy 

in electromagnetic plasma. The second contribution determining by the function S"^ jj^ 
is associated with the scattering of wave by the Debye screening shell of the particle. 
However, in contrast to the electromagnetic plasma here, the scattering is accounted for 
not through the oscillation of a coherent polarization cloud which surrounds the quark 
as a result of interaction with impingind wave, but as a consequence of the initiation of 
precession of color vectors with particles forming this cloud in the wave field uj\^. This 
process of scattering represents a pure collective effect. For calculation of its probabihty 
it is necessary to solve the kinetic equation describing a charge motion in a polarization 
cloud in the field which is equal to the sum of fields of impinding wave (|60|) and a charge 
producing a screening cloud. 



The remaining term in the right-hand side of (|6^ ) describes the interference of the 
above-mentioned scattering mechanisms. It is easily to see this having used 



_3r,2 / „,4.„ PO 



Im(-iS'fc,fcJ =TTg k2 I d p — {pu{k)){pu{ki))6{pk2){k2dpj\feq). 
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Thus, summing the preceding, instead of ( [5^ ) we have 

V(k) = - /K-OQp(k,kO(^)po^^^rfVrfki, (6^ 

J ^ Co",' ^ d»n 



where 






^uT. 



5(4 - < - v(k - ki)) I A^(k, ki) + A^(k, ki 
and the expressions for A^ and A"^ are given by ( 159|) and (0), respectively. 

Now we note that it is convenient to interpret the terms entering in the A = A^ + 
A"^ _j_ ^s ]^y ^gij2g a quantum language. In this case the term A^ connected with 
the Thomson scattering can be represented as the Compton scattering of the oscillation 
quantum (plasmon) by QGP particle. A"^*^'^^ defines the scattering of a plasmon through a 
virtual wave with the propagator l/(fc| — n'(/c2)), where a vertex of a three-wave interac- 
tion is induced by a self-action of a field. The term A'^^^^^ defines the plasmon scattering 
by a particle through a virtual wave with the same propagator and effective vertex of 
three-wave interaction connected with the medium effects. In this case A fulfils role of 
the scattering amplitude. 

The kernel (5p(k, ki) possesses two main properties. The following inequality results 
from definition (|69|) 

gp(k,ki)>o. (70) 

Next, from ( |59D it follows that A^(k, ki) = A*^(ki,k). The correctness of this equality 
results from the conservation law: cj^ — cu^;.^ — v(k — ki) = 0. For A"^ we have the 
similar relation: A'^(k, ki) = A*"^(ki, k). Its proof trivially follows from the definitions of 
Si%, Si2 ^nd e^k): 

- iSi'l = gm + kr)u{k2)){u{k)u{k,)) - {k,u{k)){u{k,)u{k2)) - 2{ku{k,)){u{k)u{k2))}^ 

(71) 

.C(^^)_ .3 [ ^4^, iP'u{k2)){p'u{k)){p'u{ki)) f {kdp>Kg) (^l^M4g)\ .„.^ 

1 + x 



e\k) = l + ^-^[l-FQ],F{x) = pn 



1 — X 



- ine{l - \x\)]. (73) 



The consequense of these equalities is a main property of a symmetry kernel Qp{k, ki) 
with respect to permutation of a wave vectors k and ki 

Qp(k,ki) = Qp(ki,k). (74) 
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Now let us consider consequence of the properties ( |70D and (0). Integrating (^) over 
(ik/(27r)^, and using (ffil), we find 

dt "dtW ^^^(2nr)-^^ 

i.e. tlie general number of plasmons A^' in the process of the nonlinear scattering conserves 
exactly. 

From ( |68l) it follows that in the case of a global equilibrium plasma, when 

dpo 

and by inequality ([70|), waves of a high frequencies are damped out, and a smaller ones are 
increased. Thus, due to the nonlinear interaction of a plasma oscillations a pumping-over 
from short to long waves occurs in the spectrum, practically leaving its total energy fixed. 
Therefore, the nonlinear decrement 7'(k) defines actually the inverse time of a spectral 
pumping. 

In the limit of |k| — i> from ( |68|) we obtain 

V(0) = - fiuj,i-ujlJQ,iOM)(^)po^^^d'pdk,. (75) 

The frequency ul^ is monotonically increasing function of a wave number, therefore 
ujpi — cj^^ < 0. From this following inequality folows 

V(o) < 0, 

i.e. |k| = 0-mode is not damped, in contrast to Ref. [8]. From the physical point of 
view this result is clear. As it was shown above, nonlinear interaction of waves leads to 
the spectral pumping from short to long waves. The mode |k| = is a limiting in this 
series, and therefore 7^(0) < 0. It is clear that in the region |k| ~ the other nonlinear 
mechanisms switching off instability are to come into effect. One of such mechanism, 
connected with regard to the terms of a higher-degree of nonlinearity in expansion of the 
color current (|T4|) will be discussed in Conclusion. 



9. THE ESTIMATE OF ^(0) 

Now we consider the quantity 7'(0). As in the paper [8] it is convenient to introduce 
the cylindrical coordinate system in which the direction of polar axis is selected to be the 
same as the direction ki. Then the coordinates for k, p, and p' are k = (|k|,a,/9), p = 
{\p\,9,ip), and p' = {\p'\,9',ip') correspondely. By and 0' we denote the corresponding 
angles between p and k, p' and k. They can be expressed as 

cos = sin 6* sin a cos{ip — /3) + cos 9 cos a. (76) 
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The expression for cos 0' is similar to above equation. 
In the hmit of Ikl — i> the kernel f|69D transforms to 



(77) 
where 

pL = «-^pO/|ki|>0. 

Using definitions (|59D and (|67|), we find the expressions A^(0,ki) and A'^(0,ki). 
Turning |k| to zero in (^9|), we have 

A2(0,ki) =/|ki|cos0K^cos^-l). (78) 

Here, we denote the phase velocity of longitudinal oscillations by f^^ = cu^yikil. 
Calculation of A'^(0, ki) is more complicated. From (|67|) we obtain 

A^(0, ki) = - -^ ^-, — — y^- — — ^-T (79) 

lkl"o Ikl ^ '^''=1 ^fc,fciJa;=wj^,wi=wJ^ • 

Here, we use definition of the function e\k) (ffSl). Dividing (|7TD by |k| and going to the 
limit |k| — >■ 0, after the simple, but slightly cumbersome computations, we define 



lim ^ "'"' 



|k|->0 Ik 



'" =2^cosa^p,|ki|3{^J-k?(l-<pLj}- (80) 



k' i k 



For definition of the limit of {—iSl ^^^)/|k|, instead of ( [7^ ) we use expression, which is 
defined from it, if the following identity is accounted for 

1 1 _ / 1 1 X 1 

p' k2 + ip^e p' ki ^p'k2 + ip'Qe p'ki' p'k 

This replacement leads to more simple limiting expression 

tS^r!'J, 3g - r2n rl 



lim ' ''''' 



|kHO Ikl 



;.a;i=w; 47r ^'7o 



UJ^UJl^. UJl=UJ, 



'li r d^' I d{cose')cos(l)'\i\{vl^ cosd'-l){pi^ cosO'-l] 



pi, ^Li 



,Pki ~ cos6'' — ie f k^ — cos 6*' ^ 
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Here, we take into account that in view of definition of the equihbrium function Meq and 
(0)(for/i = O): 



/+00 r+oo 

|pf rflp'l / p',dp\ 

-oo J —oo 



, ^j dA4,(p;)) _ 3 u^l, 

° dp^ 47r ^2 



Now we substitute instead of cos0' the expression similar ( |76D in (^TJ). Next integrating 
over Q' we obtain 



hm — — -p — 



3(7^piCosa VkiPki ^Tj^ ,/ j [Pki(l - (PkJ ) 

(l-i^(-pLj) + ^<(l -«)')]}• (82) 

Substituting (^) and (|82D in (|79|), we find unknown expression for A'^(0,ki). 
Now we return to the expression for 7^(0) (|75D, rewriting it as follows 

V(o) = |g(ki)<rfki, (83) 

where 

g(kO = ^^^^ f d^pQ,i0M)po^^^ = 144^ r^^ /' di^ose)Q,{0M). 

^ki -^ dpo 47r5f2 <^ Jo J-i 

(84) 
Let us consider the integrals over angles ip and 6 in the kernel Q(ki) (p4D. These angles 
enter into A^(0, ki) and (5-function in ([77|) , and the element A'^(0, ki) is independed from 
those at all. 

Now we explicitly select the angular dependence on square of the scattering amplitude 
entering in (5p(0, ki) 

I A^(0, ki) + A^(0, ki) p = g^ cos^ k2(t;(,^ cos 6 - 1)^+ 

+2/cos0|ki|(w(,^cose-l)ReA^(O,ki) + |A^(O,ki)|^ 
Using the formula of the angles connection (^), and taking into account relations 

p2lT p2lT 

/ cos^0(i(y9 = 7r[(3cos^6' — 1) cos^a + 1 — cos^6'] , / cos0(i(y9 = 27rcos6'cosa, 
Jo Jo 



we integrate over ip in (|8^). The remaining integral over 6 by the 5- function is computed 
elementary. Summing preceding, we find 

^C^i) = ~T~ -1~ ( ^, ./ A2 .1.2 ) ( ^:T ) , .4{(3(PkJ -l)cos«+ 
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'^i='^Li 



+ l-(pLj'}-3pl,cos2«ReA^(|ki|) + ^cos2a|A^(|ki|)r]0(l-pLj. (85) 

Here, instead of A'^(0, ki) we introduce a new function depending only on |ki| by means 
of relation 

A^(0, ki) = g^\k,\ cos« (1 - t;i^pi,JA^(|ki|). 

The pecularity of obtained expression (|85D is the absence of the angle dependence on 
f3. This enables us to represent 7'(0) in the form 

-f\0) = -2n \ki\^d\ki\ d{cosa)Q(ki)Wi^. (86) 

If we consider that the isotropy of the oscillations directions takes place at the time 
interval which is much less than the characteristic that of the nonlinear interaction, then 
the spectral density Wl__^ can be considered isotropic with respect to ki directions. This 
enables us to integrate over angle in (^) to completion. Now we introduce the spectral 
function 

W^l'k,|=47r|ki|X„ 

such that the integral /q°° Wj^rflkil = W'' is total energy of longitudinal oscillations in 
QGP. Substituting ( ^5]) into (|86D and integrating over d{cos a) we find more suitable form 
for estimate 7^(0) 

V(0) = - / Q{\k,\)Wl^^^d\k^\, (87) 

where 



2co,, [«)2-k?][3^^V(«)^-kD] 



A^dkil) 



1 



/ . U2 1 _ /,,/ \2 



< / kf 1 - {v- 



HVUi-^(-pU+^7r3fta^(i 



ki 



k? + 3^,Vi - f(-pLJ] L ^'""^^^ ' '^'''' m-ipL)')Hi pic. 



kil- 



|ki| 1-<pL/^ 1-(pU'^' ' |ki|(l-<pLj 

Here, we set 6'-function entering in (|85|) equals to unit, since the function p[^^ by its 
definition, for any |ki| satisfies the inequality 

pL < 1- 

For estimate of 7^(0) order in ( p7D we cutt-off upper integration limit on the char- 
acteristic value |ki| ^ ujpi. As was shown in [13] in the region |ki| > ujpi there is the 
strong damping of longitudinal oscillations, therefore here, this cutting has meaning. Be- 
sides, the realized above analysis of separation of leading in g terms is obeyed only in 
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the long-wavelenght spectrum region. Let us approximate the kernel (5(|ki|), using the 
approximations 

I I 3 ki I ujpi 

Leaving the leading in |ki| term in the expansion (pS]), we obtain 

Note that here, the contribution to (5(|ki|) from the Compton scattering process proves 
to be negligible by comparison with two other contributions. The function Wl__^ is approx- 
imated by its equilibrium value: Wl__^ ~ AnT and therefore 

Wfi,^^ ^ levr^lkil^T. (90) 

Substituting (^) and ( pOD in (P7|), we finally define 

V(0) ^ -m,g^T. (91) 

10. THE GAUGE DEPENDENCE 

In this section we consider the problem on the gauge dependence of nonlinear Landau 
damping rate of longitudinal oscillations. Let us compare derived expression for 7'(k) (p8| ) 
with the kernel ( |69D in covariant gauge, with similar expression computed in temporal 
gauge Aq = 0. 

As we have mentioned in Introduction, the first nonlinear Landau damping rate in 
temporal gauge was calculated in [8]. The inaccuracies in calculations were made in 
obtaining of equation for the second correction of a gauge field A^'^'{k) (formula (3.19) in 
[8]). The elimination of these inaccuracies leads to the expression 

(.;(°))^e(.;(°),k)A(^)(fc)=, ^ ^^^[A«(^i), ^^^H^]" 



-g 



d'^P Ur ^df^°Kp) rf/^°Hp)x ,, f^G'W(p) 



(27r)3EpL^^^ dEp ^ dEp )+2^^ dEp J p ■ fcW + 2/0+ ^^^^ 
E lA^ |P:kp.k,p-k, (,)(,) 

Here, we use notations accepted in [8]. The distinction ( P^ from (3.19) in [8] is as 
follows. Firstly - this is the availability of the first term in the right-hand side of (^2]), 
that is connected with self-action of a gauge field, and is of the same order in a soft region 
of excitations as the second term, connected with medium effects. As shown above this 
term also contributes to process of the nonlinear scattering of plasmons by QGP particles. 

28 



Second distinction lies in fact that instead of expression in braces in the second term in 
the right-hand side of Eq. (^) in [8] the following expression 



U)2 



p ■ kP + i/0^ 



is used. The principal point is presence of the factor 1/2 in our case. This factor enables 
us finally to lead the expression for nonlinear Landau damping rate in the temporal gauge 
to the form that is similar to (pHf), (pQ]). Further we obtain the expression for 7'(k) in 
this gauge following our reasoning, do not repeating calculations in [8]. 

In temporal gauge the Yang-Mills equation (|^) has the form 

d^F'^'^ix) - tg[A^ix), F^'^ix)] ~ r'Wu^A^ix) = -fix). 

Here, as the fixed four-vector in the gauge condition n^A'^ = we choose the four-velocity 
of plasma: n^ = u^. 



In this gauge instead of the propagator (^8]) we have 



'D^uik) = -[k'^g^,u - k^ku + i ^u^Uy - Ti^,u{k)] ^ 



^IJLV ^ av f. "' 
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In obtaining of the last equility in (|9^) we use relation [15] 



eyr^^M- (93) 



and introduce notation 



UfiUi, — -j-^Q^u -j—y—Zk u Ly^u H T2 -'-^m'^ 
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where C^^ = —{Ufj,ku + Uykfj)/y^^2k!hi?. In the rest frame of a plasma, tensor Q^^ has the 
structure 

Further we take, ^ = 0, i.e. Aq = is imposed strictly. By virtue of expansion of the 
polarization tensor II^i, = P^j/II* + Q^i^Il' and the properties of tensor structures P,Q,C 
and D [15], the expression II' remains unchanged, i.e. 

Instead of expansion of the spectral density p7|) now we have 

T — P P J^ P f) 
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Connection between the spectral density /^ and the density of the energy of longitudinal 
oscillations is unchanged. 

From the above it is easy to appreciate that for obtaining of nonlinear Landau damping 
rate in the temporal gauge it will sufficiently in the expression (7.9) to replace the projector 
Q hj Q, and in the expression for S (HH) use propagator (0^) instead of (ff3). Finally we 



obtain the same expression (|68D , where now in the kernel (|69D instead of the scattering 



amplitudes A^(k, ki), A*^ (k, ki) and A"^ (k, ki) it is necessary to use expressions 

g' (kv)(kiv) [(4)'-k2][«)'-k?] 



Ao^(k,ki 



|k||ki| c<j[. — (kv) ujI-uj- 

1 -^S^''''^^ 



Af"\kM)- ' i ' -^-okM 



|k||ki|ki \kl e\k2 



U)=Uj\^,U)l=U)\^ 



respectively. Here, 



'"-'OkM 



_^giii) ^ _ z k'klkl r, , (p-k)(p%)(p%) akd,.K,) _ {k,d,.K,) \ 
^'''''^ ^ UJUJ1UJ2 J ^ p'k2 + ip'oe V P'k P'^i J 

Now we compare the scattering amplitude A in the covariant gauge with appropriate 
scattering amplitude Aq in Ao-gauge. Using the above-mentioned explicit expressions for 
amplitudes, as the result of simple calculations we obtain 



g' ekj - u^ ul 
|k||ki| 2ujuji 



A - Ao = „,Mui o^.^. ' (k + ki)(v - v||) ^ 0, (94) 



where vy = k2(k2v)/k2. Here the term with vy in the right-hand side depends on the 
contribution to the scattering amplitudes the terms with longitudinal virtual oscillation. 
The last expression explicitly demonstrates the gauge-noninvariant character of obtained 
decrement of nonlinear Landau damping. 

11. CONCLUSION 

In our paper it was shown that the nonlinear interaction of longitudinal eigenwaves 
leads to effective pumping of energy across the spectrum sideways of small wave numbers. 
Consequence of this fact is the inequality: 7^(0) < 0, i.e. k = - mode is increased. The 
real nonlinear absorption of the energy of plasma waves by particles in QGP is the effect 
of higher-order. 

It is clear that growth of k = - mode is consequence of chosen approximation. In the 
region of a small |k| effects, described by nonlinear terms in the expansion of the color 
current of higher-order in the field, come into play . One of such possible nonlinear effects 
is known from the theory of electromagnetic plasma [11]. This is as follows. 

30 



By the effect of pumping, all plasmons will be tend to concentrate near a small |k| = 
|ko| -^ 0. However, phase space, which the plasmons are occupied, proportionaled to 
|kop will be also highly small. By virtue of this fact the intensive collision of plasmons is 
arised. This is the process 

l + h^l' + l[, 

which is to lead to the scattering of plasmons from region of a small |k| and thus to 
suppression of increase of k = 0-mode. The probability of this four-plasmon interaction 
is defined by the preceding method from the nonlinear current of fourth order jj*^^-* with 
regard to the process of interaction iteration of higher-order in field. 

Let us consider in more detail approximations scheme, which we use in this paper. In 
fact, here two various levels of the approximations are used. The first of them is connected 
with the employment of usual approach, developed in EMP to QGP, i.e. the standard 
approximation of the current in terms of the oscillations amplitude and computation of 
interacting field in the form of a series of a perturbation theory in a free field. However, 
in contrast to EMP, in our case even if the first two nonlinear orders of the color current 
are taken into account, vastly more terms, defining the nonlinear scattering of waves are 
derived. Here we use second approximation level, connected with the notions going from 
the papers on hard thermal loops or more precisely, the set of the orders estimates of a 
various terms, developed by Blaizot and lancu [5]. This set of estimates enables us to 
extract the leading terms in the coupling constant from the set of obtained terms. The 
surprising thing is that although this terms are purely non-Abelian, all basic conclusions, 
performed on the basis of these terms, coincide with the appropriate results in nonlinear 
theory of EMP in a qualitative sense. 



The obtained expression of nonlinear Landau damping rate (|68D, (pOf ) is not gauge 
invariant, that it was explicitly shown in Sec. 10. Such gauge dependence of the nonlinear 
damping rate is a specific character of the non-Abelian theory. In the Abelian plasma 
a similar problem does not arise, since the theory of the nonlinear processes in EMP 
is stated in terms of the gauge-invariant electric and magnetic fields only. However in 
the case of QGP it is impossible to avoid the work with potentials, since ones explicitly 
appear in the kinetic equations (^) and YM Eq. (^). Therefore the initial function in 
our consideration is the two-point correlative function of the form ( p3D (more precisely, 
its longitudinal part, since in this paper we have restricted ourselves to the longitudinal 
excitations in QGP only), weakly inhomogeneous and weakly nonstationary, having no 
directly physical meaning and being the gauge-noncovariant value. 

By virtue of fact that the process of the nonlinear scattering of longitudinal waves is 
the physical process, the scattering amplitude A, entering in the kernel (|^) is bound to be 
gauge-invariant (within used in this paper approximation the spectral density of energy 
Wl_^ is gauge-invariant value). However, as it was shown in Sec. 10 (Eq. (|94D ), the scattering 
amplitudes calculated in the covariant and Ag-gauges are not coincident. The simplest 
analysis of the right-hand side of (p^) points to the fact that difference A — Aq is equal 
to zero if in the last parentheses in the right-hand side besides longitudinal in relation 
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to the vector k2 component of velocity vy the transverse component v_l: (v_Lk2) = is 
present. Intuitively clear that such term is arised from discarded contribution to nonlinear 
scattering of term with transverse virtual oscillation. Howerever if the longitudinal and 
transverse virtual oscillations taken into account, already it is impossible to write 7'(k) 
in compact and transparent form identical with those of (|68D , (^9D . Therefore in this case 
checking of gauge invariance becomes a very nontrivial problem and it is the subject of 
specific research. 



Acknowledgements 

This work was supported by the Russian Foundation for Basic Research (Project No.97- 
02-16065). 



32 



References 

[I] E. Braaten and R.D. Pisarski, Nucl. Phys. B337, 569 (1990); J. Frenkel and J.C. 
Taylor, Nucl. Phys. B334, 199 (1990). 

[2] U. Heinz, Phys. Rev. Lett. 51, 351 (1983); Ann. Phys. (N.Y.) 161, 48 (1985); 168, 
148 (1986); H.-Th. Elze, M. Gyulassy, and D. Vasak, Phys. Lett. B177, 402 (1986); 
Nucl. Phys. B276, 706 (1986); H.-Th. Elze and U. Heinz, Phys. Rep. 183, 81 (1989); 
St. Mrowczynski, Phys. Rev. D39, 1940 (1989); Yu. Markov and M. Markova, Teor. 
Math. Fiz. Ill, 263 (1997) (Theor. Math. Phys. Ill, 601 (1997)). 

[3] V.P. Silin, Zh. Eksp. Teor. Fiz. 38, 1577 (1960) [Sov. Phys. JETP 11, 1136 (I960)]. 

[4] P.F. Kelly, Q. Liu, C. Lucchesi, and C. Manuel, Phys. Rev. D50, 4209 (1994). 

[5] J.P. Blaizot and E. lancu, Phys. Rev. Lett. 70, 3376 (1993); Nucl. Phys. B417, 608 
(1994). 

[6] J.P. Blaizot and E. lancu, Phys. Rev. Lett. 72, 3317 (1994); Phys. Lett. B326, 138 
(1994). 

[7] E. Braaten and R.D. Pisarski, Phys. Rev. D42, 2156 (1990). 

[8] Z. Xiaofei and L. Jiarong, Phys. Rev. C52, 964 (1995). 

[9] U. Heinz and P.J. Siemens, Phys. Lett. B158, 11 (1985). 

[10] A.G. Sitenko, Fluctuations and Non-linear Wave Interactions in Plasma (Naukova 
Dumka, Kiev, 1977) [English transl. publ. by Pergamon, Oxford, 1990]. 

[II] B.B. Kadomtsev, Plasma Turbulence (Academic Press, New York, 1965); V.V. Pusto- 
valov and V.P. Silin, Proc. P.H. Lebedev Inst. 61, 42 (1972); V.N. Tsytovich, Non-linear 
Effects in Plasma (Nauka, Moscow, 1967) [English transl. publ. by Plenum, Oxford, 
1970]; V.N. Tsytovich, Theory of Turbulent Plasma (Nauka, Moscow, 1971) [English 
transl. publ. by Plenum, Oxford, 1977]. 

[12] H.-Th. Elze, Z. Phys. C38, 211 (1988). 

[13] O.K. Kalashnikov and V.V. Khmov, Yad. Fiz. 31, 1357 (1980) [Sov. J. Nucl. Phys. 
31, 699 (1980)]; V.V. Klimov, Zh. Eksp. Teor. Fiz. 82, 336 (1982) [Sov. Phys. JETP 
55, 199 (1982)]. 

[14] H.A. Weldon, Phys. Rev. D26, 1394 (1982). 

[15] K. Kajantie and J. Kapusta, Ann. Phys. (N.Y.) 160, 477 (1985); U. Heinz, K. Ka- 
jantie, and T. Toimela, Ann. Phys. (N.Y.) 176, 218 (1987). 

[16] R. Efraty and V.P. Nair, Phys. Rev. D47, 5601 (1993); R. Jackiw and V.P. Nair, 
Phys. Rev. D48, 4991 (1993). 



33 



